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Abstract. In this manuscript we provide a family of lower bounds on the indirect Coulomb 
energy for atomic and molecular systems in two dimensions in terms of a functional of the 
single particle density with gradient correction terms. 

cV 
o 

pH ■ 1. Introduction 

Since the advent of quantum mechanics, the impossibility of solving exactly problems in- 
t^. ! volving many particles has been clear. These problems are of interest in such areas as atomic 
and molecular physics, condensed matter physics, and nuclear physics. It was, therefore, nec- 
' essary from the early beginnings to estimate various energy terms of a system of electrons 
as functionals of the single particle density p^(x), rather than as functionals of their wave 
J j function ip. The first estimates of this type were obtained by Thomas and Fermi in 1927 
(see [jnj] for a review), and by now they have given rise to a whole discipline under the name 
of Density Functional Theory (see, e.g., M and references therein). In Quantum Mechanics 
of many particle systems the main object of interest is the wavefunction ip £ /\ L 2 (IR 3 ), 
(the antisymmetric tensor product of L 2 (IR 3 )). More explicitly, for a system of N fermions, 
ip(xi, ... ,Xi, ....... Xj, .. . xn) = —ip(xi, . . . , Xj, ....... Xi, ... Xjsr), in view of Pauli's Exclu- 



> 



OA ' sion Principle, and J RN \ip\ 2 dx\ . . . dx n = 1. Here, x« £ M 3 denote the coordinates of the i-th 
particle. Fi 
density) as 



particle. From the wavefunction ip one can define the one-particle density (single particle 



in 

(N ■ Pi>{ x ) = N I \ip{x,x 2 ,...,x N )\' 2 dx 2 ...dx N , (1) 

VlR3(iV-l) 

and from here it follows that L 3 p^(x) dx = N, the number of particles, and p^(x) is the 
density of particles at x £ M 3 . Notice that since ip is antisymmetric, \ip\ 2 is symmetric, and it 
is immaterial which variable is set equal to x in ([!]). 

In Atomic and Molecular Physics, given that the expectation value of the Coulomb attrac- 
tion of the electrons by the nuclei can be expressed in closed form in terms of Pip(x), the 
interest focuses on estimating the expectation value of the kinetic energy of the system of 
electrons and the expectation value of the Coulomb repulsion between the electrons. Here, we 
will be concerned with the latest. The most natural approximation to the expectation value 
of the Coulomb repulsion between the electrons is given by 

If 1 

D (p^,Pi>) = ^ / pA x ) \ _ \ Pi>{y)& x &yi (2) 
^ j \ x y\ 

which is usually called the direct term. The remainder, i.e., the difference between the expec- 
tation value of the electronic repulsion and D(p^,p^), say E, is called the indirect term. In 

1930, Dirac |J gave the first approximation to the indirect Coulomb energy in terms of the 

l 
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single particle density. Using an argument with plane waves, he approximated E by 

E ~ -c D [ pf dx, (3) 



where c D = (3/4)(3/tt) 1 / 3 « 0.7386 (see, e.g., |0|, p. 299). Here we use units in which the 



absolute value of the charge of the electron is one. The first rigorous lower bound for E was 
obtained by E.H. Lieb in 1979 using the Hardy-Littlewood Maximal Function |27|. There 
he found that, E > —8.52 J pf dx. The constant 8.52 was substantially improved by E.H. 
Lieb and S. Oxford in 1981 [PJ, who proved the bound 



E>-C J pf dx, (4) 
with C = Clq = 1.68. In their proof, Lieb and Oxford used Onsager's electrostatic inequality 



22| , and a localization argument. The best value for C is unknown, but Lieb and Oxford 
| 16| proved that it is larger or equal than 1.234. The Lieb-Oxford value was later improved 
to 1.636 by Chan and Handy, in 1999 f|. Since the work of Lieb and Oxford ||16[| , there 
has been a special interest in quantum chemistry in constructing corrections to the Lieb- 
Oxford term involving the gradient of the single particle density. This interest arises with 
the expectation that states with a relatively small kinetic energy have a smaller indirect part 

see, e.g., jy], 24, ^HJ and references therein). Recently, Benguria, Bley, and Loss obtained 



an alternative to (|4j), which has a lower constant (close to 1.45) to the expense of adding a 
gradient term (see Theorem 1.1 in ||), which we state below in a slightly modified way, 

Theorem 1.1 (Benguria, Bley, Loss 0). For any normalized wave function ip(xx, . . . , xn) 
and any 0<e< 1/2 we have the estimate 



E^)> -1.4508 (1 + e) / pfdx - \p\^) 



(5) 



where 

(VPAP\V~P) ■■= [ \?P(k)\ 2 \^k\dk = ^ [ [ \^-^ )? dxdy. (6) 
iR3 2tH J R 3 V \x - y\ A 

Here, f(k) denotes the Fourier-transform 

f(k) = f e- 2mk x f(x)dx . 



Remarks, i) For many physical states the contribution of the last two terms in is small 
compared with the contribution of the first term. See, e.g., the Appendix in 0; 

ii) For the second equality in ^) see, e.g., fL5| , Section 7.12, equation (4), p. 184; 

Hi) It was already noticed by Lieb and Oxford (see the remark after equation (26), p. 261 
on that somehow for uniform densities the Lieb-Oxford constant should be 1.45 instead 

0/1.68; 

iv) In the same vein, J. P. Perdew |23|, by employing results for a uniform electron gas in 



its low density limit, showed that in the Lieb-Oxford bound one ought to have C > 1.43 (see 
also 0;. 



After the work of Lieb and Oxford |T6[ many people have considered bounds on the indirect 
Coulomb energy in lower dimensions (in particular see, e.g., |10] for the one-dimensional case; 



21], [25], and |26] for the two-dimensional case, which is important for the study of 
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quantum dots). Recently, Benguria, Gallegos, and Tusek gave an alternative to the Lieb- 
Solovej-Yngvason bound JTB), with a constant much closer to the numerical values proposed 
in OH (see also the references therein) to the expense of adding a gradient term: 



Theorem 1.2 (Estimate on the indirect Coulomb energy for two dimensional atoms ||). 
Let ip G L 2 (M. 2N ) be normalized to one and symmetric (or antisymmetric) in all its variables. 
Define 



p^(x) = N I (x, x 2 , ■ ■ ■ , x N ) dx 2 ... dx 



AT- 



(N-l) 



If p^ G L 3 / 2 (R 2 ) and |VpJ, /4 | G L 2 {R 2 ), then, for alle>0, 

N 



E^) = (V,X)|a:i-x i |-V>-^(p^/v)>-(l + e)/8 / pf dx ~ 7T f \ V pf? dx ( 7 ) 



with 



/4\ 3/2 . 

/3 = f - J V57T - 1 ~ 5.9045. 



Remarks. TTie constant (3 ~ 5.9045 m (W is substantially lower than the constant Clsy 
481.27 found in [IS (see equation (5.24) of lemma 5.3 in 



iij Moreover, the constant (3 is close to the numerical values (i.e., ~ 1.95 ^ of [^5] (and refer- 
ences therein), but is not sharp. 

In the literature there are, so far, three approaches to prove lower bounds on the exchange 
energy, namely: 

i) The approach introduced by E.H. Lieb in 1979 [|13[], which uses as the main tool the Hardy 



Littlewood Maximal Function |27]. This method was used in the first bound of Lieb ||T3| . 



Later it was used in to obtain a lower bound on the exchange energy of two-dimensional 
Coulomb systems. It has the advantage that it may be applied in a wide class of problems, 
but it does not yield sharp constants. 



ii) The use of Onsager's electrostatic inequality []22| together with localization techniques 



introduced by Lieb and Oxford |fLq| . This method yields very sharp constants. It was used 
recently in [0 to get a new type of bounds including gradient terms (for three dimensional 
Coulomb systems). In some sense the constant 1.4508 recently obtained in is the best 
possible (see the comments after Theorem |1 . 1| ) . The only disadvantage of this approach is 
that it depends on the use of Onsager's electrostatic inequality (which in turn relies on the fact 
that the Coulomb potential is the fundamental solution of the Laplacian). Because of this, it 
cannot be used in the case of two-dimensional atoms, because l/\x\ is not the fundamental 
solution of the two-dimensional Laplacian. 

iii) The use of the stability of matter of an auxiliary many particle system. This idea was 
used by Lieb and Thirring |19| to obtain lower bounds on the kinetic energy of a systems of 
electrons in terms of the single particle density. In connection with the problem of getting 
lower bounds on the exchange energy it was used for the first time in |3J, to get a lower bound 
on the exchange energy of two-dimensional Coulomb systems including gradient terms. This 
method provides very good, although not sharp, constants. 
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As we mentioned above, during the last twenty years there has been a special interest in 
quantum chemistry in constructing corrections to the Lieb-Oxford term involving the gradi- 
ents of the single particle density. This interest arises with the expectation that states with a 
relatively small kinetic energy have a smaller indirect part (see, e.g., fTT , |24| , 28] and references 



therein). While the form of leading term (i.e., the dependence as an integral of p 4//3 in three 
dimensions or as an integral of p 3//2 in two dimensions) is dictated by Dirac's argument (using 
plane waves), there is no such a clear argument, nor a common agreement concerning the 
structure of the gradient corrections. The reason we introduced the particular gradient term, 
J R 2 | V/^/ 4 | 2 dx in our earlier work ||], was basically due to the fact that we already knew the 
stability of matter arguments for the auxiliary system. However, there is a whole one parame- 
ter family of such gradient terms that can be dealt in the same manner. In this manuscript we 
obtain lower bounds including as gradient terms this one-parameter family. One interesting 
feature of our bounds is that the constant (5 in front of the leading term remains the same (i.e., 
its value is independent of the parameter that labels the different possible gradient terms), 
while the constant in front of the gradient term is parameter dependent. 

Our main result is the following theorem. 

Theorem 1.3 (Estimate on the indirect Coulomb energy for two dimensional atoms). Let 

1 < 7 < 3 ; and a = (3 - 7)/(2 7 ). Assume, p^ G L 3 / 2 (R 2 ) and |Vp$| G L 7 (M 2 ). Let 
C(p) = 2 1 ~ p / 2 , for < p < 2 while C(p) = I, for p > 2. Then, for all e > we have, 



n r r 

E(^) = (^,V|x J -x,|-V)-^(p^P^)>-& 2 / pfdx-a 2 \Vp^dx. 
7^ J^ 2 J^ 2 



(9) 



Here, 



4X 3/2 



b 2 =[-) v / 5^T(l + e ) = /3(l + e) (10) 



where (3 is the same constant that appears in ^). Also, 



. 2 _2^C( 7 ) fl>y-l„f 7 



In particular, we have (with a fixed e) 

^2 

I7-M+ 



a 2 U^ 1+ = a/2. 



Remarks, i) Our previous Theorem is a particular case of Theorem for the value 
7 = 2, a = 1/4. 

ii) Notice that b 2 is independent of'j, and it is therefore the same as in 

Hi) The constant in front of the gradient term depends on the power 7 and, of course, on e. 
However, as 7 — > 1+, this constant converges to s/2 independently of the value of e. 



In the rest of the manuscript we give a sketch of the proof of this theorem, which follows 
closely the proof of the particular result [L2| in [|J . 
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2. Auxiliary lemmas 

First we need a standard convexity result. 

Lemma 2.1. Let x,y G R ; and p > 0. Then 

\x\ p + \y\ p < C(p)\x + iy\ p , 

where C(p) = 2 x ~ p l 2 for < p < 2, and C(p) = 1 for p > 2. The constant C(p) is sharp. 

Proof. If p > 2, the assertion follows, e.g., from the fact that Z p -norm is decreasing in p. On 
the other hand, for < p < 2, the assertion follows from the concavity of the mapping t — > t 1 ' n 
for t > and n > 1. □ 

The next lemma is a generalization of the analogous result introduced in and used in 
the proof of Theorem \1.2\ above (see 0). This lemma is later needed to prove a Coulomb 
Uncertainty Principle. 

Lemma 2.2. Let Dr stands for the disk of radius R and origin (0, 0). Moreover letu = u(\x\) 
be a smooth function such that u(R) = and 1 < 7 < 3. Then the following uncertainty 
principle holds 



Dr- 



[2u(\x\) + \x\u\\x\)]f(x) 1,a 

1/7 



< 



<\{c{l)J D \Vf(x^dx^j J (c(S)J d \x\ 5 \u(\x\)\ 5 \f(x)\ 3 ^dx 



1/6 



where 



1 27 

a 3 — 7 ' 

Proof. Set gj(x) = u(\x\)xj. Then we have, 



1 1 

- + x = 1 - 

7 



(12) 



[2u(\x\) + \x\u'{\x\)]f{x) x ' a dx = J2 [d jgj {x)]f{x) lla dx 



j=1 JD R 



1 

In the last equality we integrated by parts and made use of the fact that u vanishes on the 
boundary 8Dr. Next, the Holder inequality implies 



d r 

1 

a \Jd 



[2u(\x\) + \x\u'(\x\)]f(x) 1/a 



< 



2 \ l h 

Y,\d 3 f{x)Vdx 

R j = l 



1/5 



D 



J2\xj\ 5 H\x\)\ S \f{x)\ (1/a - 1)S dx 



R j = l 



The rest follows from Lemma 2.1 



□ 
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3. A STABILITY RESULT FOR AN AUXILIARY TWO-DIMENSIONAL MOLECULAR SYSTEM 

Here we follow the method introduced in [|J. That is, in order to prove our Lieb-Oxford 
type bound (with gradient corrections) in two dimensions we use a stability of matter type 
result on an auxiliary molecular system. This molecular system is an extension of the one 
studied in J3J, which was adapted from the similar result in three dimensions discussed in H 
(this last one corresponds to the zero mass limit of the model introduced in |7|, |?| ||). We begin 
with a typical Coulomb Uncertainty Principle which uses the kinetic energy of the electrons 
in a ball to bound the Coulomb singularities. 

Theorem 3.1. For every smooth non-negative function p on the closed disk Dr C M 2 , and 
for any a, b > we have 

aba [ ( 1 L-^] P (x)dx f | V p(*)Td* + ^ f p^dx, 

Jd r \\x\ RJ 7 Jd r o J Dr 

where 1 < 7 < 3, and a and 5 are as in fln\). 

Proof. In Lemma |2~2] we set u(r) = 1/r — l/R and / = p a . The assertion of the theorem then 
follows from Young inequality with coefficients 7 and 5. □ 

And now we introduce the auxiliary molecular system through the "energy functional" 



i{p) = a 2 [ \Vp a \^dx + b 2 [ p 3/2 dx- [ V(x)p(x)dx + D(p,p) + U, 
Jr 2 Jr 2 Jr 2 



(13) 



where 

K 



V ^ = f2r^lTr D (P>P) = \f p(x)-^—p(y)dxdy, u= 

T^\ x -Ri\ 2i t 2 xR2 \x-y\ ^f-i> 



z 2 



I Ri — Rj I 



i=l 1 11 1 31 l<i<j<K 

with z > and Ri G M 2 . As above we assume 1 < 7 < 3, and a = (3 — 7)/ (27). The choice 
of a (in terms of 7) is made in such a way that the first two terms in ( |T3"D scale as one over a 
length. Indeed, let us denote 



K(p) = a 2 [ |VpT dx + b 2 [ p 3 / 2 dar. 
Jr 2 Jr 2 



Given any trial function p e L 1 (M 2 ) and setting p\(x) = X 2 p(Xx) (thus preserving the L 1 
norm), it is simple to see that with our choice of a we have K(p\) = XK(p). 
If we now introduce constants a, b±, b 2 > so that 



~ 2 _ a 7 C( 7 ) 



a* = w/ (14) 
2«7 v ; 

b =^aJT + b * 

(again with 5 given by fljjD ), we may use the proof of [U Lemma 2.5] step by step. In 
particular, 

£(p) > b\ [ p 3 / 2 dx - [ Vp dx + ab 2 I ( 9| p I " pWx + D(p, p) + U, 
Jr 2 Jr 2 ~[Jb 3 \ 2 \ x - R j I D j J 

where 

Dj = - mm{\R k -R j \\k^ j}, 
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and Bj is a disk with center Rj and of radius Dj. 
Thus as in ||, Lemma 2.5] we have that, for 



it holds 



K 



3=1 1 



z < ab 2 /2 
z 2 4 



(2z 3 (n - 1) + 7ra 3 & 3 ) 



(15) 
(16) 



8 276? 

Consequently we arrive at the following theorem. 

Theorem 3.2. For all non-negative functions p such that p E L 3 / 2 (M 2 ) and |Vp a | G L 7 (M 2 ) ; 
we have that 

£0) > o, (17) 

provided that 

(18) 



z < max h(a) 
<re(o,l) 



. a I ~ 2 3 - 7 
^((jj = mm < — | b 



C 



7 



( 7 -l)/ 7 



'I -a) 



27 6 4 



(19) 



2 \ " 7- 1 ~ V 7 - iy v " ' / ' 64 5vr - 1 

wit/i a g'ifen by ^L^)- 

In order to arrive at (|l^) we set 62 in flTBD to be the smallest possible under the condition 



(|15|), i.e., 62 = 2z/a, and we introduced a = byb . 



4. Proof of Theorem 1.3 



In this Section we give the proof of the main result of this paper, namely Theorem |L3[ We 
use an idea introduced by Lieb and Thirring in 1975 in their proof of the stability of matter 



|T9| (see also the review article |T2| and the recent monograph |T7|]). This idea was first used 
in this context in [Q • 

Proof of Theorem \1.3j . Consider the inequality (|i~7D , with K = N (where N is the number of 
electrons in our original system), z = 1 (i.e., the charge of the electrons), and Ri = Xi (for all 
i — 1, . . . ,N). With this choice, according to (|18"D, the inequality ([[7]) is valid as long as a 
and b (that are now free parameters) satisfy the constraint, 

1 < max h(a) (20) 

0-6(0,1) 

with cr (which maximizes h(a)) such that h(a ) = 1. Let us introduce e > and set 



do = 1/(1 + e)- Then the smallest b such that the assumptions of Theorem |3.2| may be in 
principle fulfilled reads 



V57r-l(l + e) 



Hence a has to be chosen large enough, namely such that 



3/2 



(21) 



- ( b 2 3 ~ 7 



C 



7 — 1 V 7 — 1 



7 



(7-l)/7 



8 



BENGURIA AND TUSEK 



which due to ([T|) implies 

3/2 ^ \ T" 1 



a 2 



Since 



lim C( 7 ) = v 7 ^, lim — ic(—^— X) = 1, 



we have (with a fixed e) 



a 2 | 7 ^i+ = \/2. 



Then take any normalized wavefunction if)(xi, £2, • • • , xn), and multiply (|17]) by \ip{xi, ■ ■ ■ > x n) | 2 
and integrate over all the electronic configurations, i.e., on ¥L 2N . Moreover, take p = p^(x). 
We get at once, 

E{$) = (^y^lxi-Xjl-^-D^p^,) > -a 2 [ \V P ydx-b 2 [ p^dx (23) 

~^ Jr 2 Jr 2 

provided a and b satisfy ( p2|) and (|2l|) , respectively. □ 

Remark 4.1. In general the two integral terms in (Qj are not comparable. If one takes a 
very rugged p, normalized to N , the gradient term may be very large while the other term can 
remain small. However, if one takes a smooth p, the gradient term can be very small as we 
illustrate in the example below. Let us denote 

L(p) = [ p(x) 3/2 dx 
Jr 2 

and 

G{ P )= I (\v P (xr\ydx. 

Js 2 

with a = (3 — 7)/ (27) . We will evaluate them for the normal distribution 

p(\x\) = Ce- A|x|2 

where C, A > 0. Some straightforward integration yields 

2tt 



L = C Z ' 2 - 

3 A' 



while, 

With C = NA/tt, 
and we have 



G = C a7 7r2 7 (Aa) (7/2) - 1 r (l + Tj 7-Ct/ 2 )- 1 . 
p(\x\) dx = N, 



L 

i.e., in the "large number of particles" limit, the G term becomes negligible, for all 1 < 7 < 3. 
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